Multiplicity distribution and source deformation in full-overlap U+U collisions 



Anthony Kuhlman and Ulrich Heinqj 
Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

(Dated: February 9, 2008) 

We present a full Monte Carlo simulation of the multiplicity and eccentricity distributions in U+U 
collisions at sfs = 200 A GeV. While unavoidable trigger inefficiencies in selecting full-overlap U+U 
collisions cause significant modifications of the multiplicity distribution shown in 0, a selection of 
source eccentricities by cutting the multiplicity distribution is still possible. 

PACS numbers: 25.75.-q, 25.75.Nq, 12.38.Mh, 12.38.Qk 
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In we advocated the use of full-overlap collisions be- 
tween deformed uranium nuclei to probe open questions 
at the Relativistic Heavy Ion Collider (RHIC). Specif- 
ically, such collisions can be used to test the hydrody- 
namic behavior of elliptic flow to much higher energy den- 
sities than currently possible with non-central Au+Au 
collisions, and the large and strongly deformed reaction 
zones produced in such collisions will allow for a detailed 
examination of the path length dependence of the energy 
lost by a fast parton as it travels through the plasma 
created in the collision. 

The calculations presented in were based on the 
assumption that full-overlap U+U collisions can be effi- 
ciently triggered on by using the zero degree calorimeters 
(ZDCs) of the RHIC experiments to discriminate against 
collisions with spectator nucleons flying down the beam 
pipe. By cutting the multiplicity distribution for the 
thus selected full-overlap collisions one can further se- 
lect subevent classes with different spatial deformations 
of the created fireballs. In this short note we explore 
these assumptions in more quantitative detail, using a 
Monte Carlo simulation of the distributions of specta- 
tor nucleons and charged particle multiplicity for U+U 
collisions for arbitrary impact parameter and relative ori- 
entation between the two deformed uranium nuclei. 

Our calculations are based on a Glauber model 
parametrization of the initial entropy production in these 
collisions, with standard Woods-Saxon form for the den- 
sity distributions of the colliding nuclei (see 0] for de- 
tails). [Possible modifications arising from Color Glass 
Condensate (CGC) jjlinitial conditions Q will be shortly 
discussed at the end The initial entropy density in 
the transverse plane at z = (with z denoting the beam 
direction) is determined by a combination of terms pro- 
portional to the wounded nucleon (n wn ) and binary col- 
lision (ribc) densities: 

s(r_L;$) = K s [«n wn (r i ;$) + (h«)n bc (r i ;$)]. (1) 
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FIG. 1: Multiplicity per participant pair versus the number of 
participants Apart = J d 2 rj_n vn (rj_) for 200AGeV Au+Au 
and U+U collisions. Solid line: Glauber model fit to the 
Au+Au data from PHOBOS 0. Dashed line: prediction for 
full-overlap U+U collisions. 



is adjusted to reproduce the charged particle multiplic- 
ity density dN^/dy measured at midrapidity in central 
200 A GeV Au+Au collisions at RHIC |5j , assuming pro- 
portionality of dNch/dy with the total entropy produced 
in the transverse plane. The wounded nucleon scaling 
fraction is tuned to a = 0.75 to reproduce the central- 
ity dependence of dN^/dy (see Figure After fitting 
n s and a to the Au+Au data, we use the same parame- 
ters to predict the multiplicities for U+U collisions at the 
same y/s. The results for full-overlap U+U collisions are 
shown in Figure^by the dashed line; low (high) multipli- 
cities correspond to side-on-side (tip-on-tip) collisions as 
indicated, due to their smaller (larger) binary collision 
contribution. 

To determine the multiplicity distribution, we intro- 
duce Gaussian event-by-event fluctuations of the multi- 
plicity n = dNch/dy via Q 



Here $ is the angle between the beam direction and sym- 
metry axis of one of the two U nuclei; in full-overlap col- 
lisions the symmetry axis of the other U nucleus lies in 
the same plane and forms an angle of <& or 7r — <& with 
the beam axis. The normalization constant k s in Q 
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where n($) is the average charged particle multiplicity 
computed from Eq. in a U+U collision with orien- 
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tation angle $, and a width of a = 0.6 has been shown 
to yield good agreement with PHOBOS data Q. The 
multiplicity distribution is then obtained by integrating 
over 4>. The resulting distribution [l| is shown by the 
gray line in Figure below. Its double hump structure 
results from the Jacobian dn/d$>, and its asymmetry is a 
consequence of the fluctuation width being proportional 
to the mean multiplicity n. Note that the non-linear de- 
pendence of the charged multiplicity on the number of 
participant (wounded) nucleons, arising from the binary 
collision component in our parametrization |JT)|. leads to 
a ~ 15% variation of the charged particle multiplicity 
among full-overlap U+U collisions as the relative orien- 
tation of their symmetry axes is varied over the accessible 
range. 

The calculations of multiplicity and eccentricity dis- 
tributions presented in Q] rely on the assumption that, 
by monitoring spectator neutrons in the backward and 
forward zero degree calorimeters (ZDCs) of the RHIC 
experiments, full-overlap collisions can be perfectly dis- 
tinguished from those collisions where the two nuclei are 
slightly misaligned. This is impossible in practice since 
even fully aligned collisions in general have a small num- 
ber of spectator nucleons, arising from the dilute nuclear 
surface, and this number is larger for side-on-side than 
for tip-on-tip collisions (see Figure^). Therefore, slightly 
misaligned tip-on-tip and fully aligned side-on-side colli- 
sions can have the same -/V part and the same ZDC signal. 
To assess the contamination from collisions with imper- 
fect overlap on the multiplicity and eccentricity distri- 
butions requires a more comprehensive study which in- 
cludes non-central U+U collisions. This is the point of 
this short note. 

A general U+U collision is parametrized by 5 pa- 
rameters, the impact parameter b and two Euler angles 
fl = ( < E > ,/3) for each nucleus describing the orientation of 
its symmetry axis relative to the beam axis and impact 
parameter direction. Equation Q for the initial entropy 
density must therefore be generalized to 

s(r ± ; 6,n 1 ,f2 2 ) = k s [an wn (r±; b, Qi, Q 2 ) 

+ (l-a)n bc (r±;6,ni,n 2 )]. (3) 

The multiplicity distribution is then calculated from 
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Evaluating this 5-dimensional integral by Monte Carlo 
integration, we obtain the multiplicity distribution shown 
in Figure [21 Its right-most part contains the full-overlap 
collisions. 

We can now try to select the latter from the overall 
event population by placing stringent cuts on the num- 
ber of spectators (= 2 x 238 — N palt )- The distribution 
of the number of spectators is shown in Figure below. 
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FIG. 2: Multiplicity distribution (normalized to unit total 
probability) for x /s = 200AGeV U+U collisions. The distri- 
bution was generated from an impact parameter weighted sum 
of approximately 200,000 events with all possible orientations 
and impact parameters between the two nuclei. 



In Figure|21we show the multiplicity distributions associ- 
ated with the 5% and 0.5%, respectively, of events with 
the lowest spectator counts 0- It is immediately ob- 
vious that contamination from slightly misaligned colli- 
sions is sufficient to completely destroy the double-hump 
structure of the ideal full-overlap case, replacing it with 
a single peak. By selecting low-spectator events, we bias 
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FIG. 3: Multiplicity distribution for the ideal case of only full- 
overlap collisions (gray line), and for the 0.5% and 5% U+U 
collision events with the lowest spectator counts (solid blue 
and dotted red histograms, respectively ). All distributions 
are normalized to unit total probability. The vertical lines 
indicate 5% and 95% cuts on the blue multiplicity distribution 
for the 0.5% spectators event class. 
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the symmetry axes of the two nuclei approximately par- 
allel. This suppresses the contribution from side-on-side 
configurations under the left peak of the idealized double- 
hump structure. At the same time, slightly misaligned 
tip-on-tip collisions fill in the dip between the two humps 
from the idealized case 0. The result is a single-peaked 
multiplicity distribution whose center moves left (towards 
lower multiplicities) as the cut on the number of specta- 
tor nucleons is loosened. 

Nevertheless, for sufficiently tight spectator cuts, we 
still expect the collision events corresponding to the left 
edge of the multiplicity distributions shown in Figure0]to 
have a larger contribution from strongly deformed side- 
on-side collisions than the events from the right edge 
(which will be mostly tip-on-tip collisions with small or 
zero source eccentricity). Following our previous sugges- 




x 



FIG. 4: Normalized eccentricity distributions, after making 
cuts on the multiplicity distributions for the two event classes 
shown in Figure|3](top panel: 5% spectator cut; bottom panel: 
0.5% spectator cut). Looser spectator cuts lead to broader 
eccentricity distributions and less well-defined average eccen- 
tricities. 

tion to select source eccentricities by cutting the multi- 
plicity distribution of "zero spectator" collisions, we per- 
form such cuts on the more realistic distributions shown 
in Figure Figure 01 shows that it still possible in this 
way to select event classes with a given average source ec- 
centricity: By taking the 0.5% of events with the lowest 
spectator count from Figure 01 (solid histogram) and cut- 
ting once more on the 5% of events with the lowest mul- 
tiplicity, we obtain the eccentricity distribution shown 
by the black histogram in the bottom panel of Figure 01 
This event class has an average source deformation e x 
of about 18%, corresponding to Au+Au collisions with 
impact parameters around 5.5 fm. On the other hand, 
taking the same 0.5% spectator cut and selecting the 5% 
events with the largest multiplicities we obtain for the ec- 



centricity distribution the gray histogram in the bottom 
Figure0] this distribution peaks at e x = and has a very 
small average spatial deformation. 

If one loosens the spectator cut to 5% instead of 0.5% 
(dotted histogram in Figure |3J) and performs the same 
multiplicity selections (5% lowest or largest multiplici- 
ties, respectively), one obtains the eccentricity distribu- 
tions shown in the top panel of Figure 01 Clearly, these 
distributions are much broader than with the tighter 
spectator cut, and the average eccentricities shift down 
from 17.7% to 14.2% for the low-multiplicity selection 
and up from 2.2% to 4.3% for the high-multiplicity selec- 
tion. Note that, since the looser spectator cut allows for 
an increased contribution from non-zero impact param- 
eters, the eccentricity of the nuclear overlap region can 
actually exceed the ~ 25% ground state deformation of 
the single-uranium density distribution projected on the 
transverse plane. This gives rise to the right tail of the 
black histogram in the top panel of Figure 01 A typical 
event from this tail is shown in Figure One sees that 
the 5% spectator cut allows for sizeable nonzero impact 
parameters and numbers of spectator nucleons, and that 
very tight ZDC cuts are required to ensure almost full 
overlap of the two uranium nuclei. 
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FIG. 5: Distribution of spectator nucleons for 200AGeV 
U+U collisions The inset shows initial entropy density 
contours (s = 20, 40, 60, 80, 100 fm -3 from the outside in) for 
a high-eccentricity event (e x = 0.325) which is included in a 
"5% spectator cut" event sample, as indicated by the arrow. 
The impact parameter of this collision is b = 3.14fm, and it 
produces N spec = 87 spectator nucleons. 

The detailed shapes of the eccentricity distributions 
shown in Figure 01 are expected to depend somewhat on 
our parametrization 11I3|I of the initial transverse den- 
sity distribution of the produced matter. It was shown 
in Ref. ,4f that initial conditions motivated by the Color 
Glass Condensate picture of low- a; gluon saturation in 
large nuclei at high energies [jj produce transverse den- 
sity distributions which fall off more steeply near the edge 
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than the more Gaussian-like distributions resulting 
from our Eqs. (11 LSI) . This might result in somewhat larger 
eccentricities and narrower eccentricity distributions (i.e. 
better defined average eccentricities) than those shown in 
Figure 0] However, significantly higher statistics would 
likely be needed to clearly see such differences. 

These results show that the suggestions made in Q 
for using U+U collisions to explore in more detail the 
ideal fluid dynamic nature of elliptic flow and the path 
length dependence of radiative parton energy loss are rea- 
sonably robust against trigger inefficiencies, and that a 
meaningful U+U collision program at RHIC is, in fact, 
feasible |l2^ . Simultaneous strict cuts on small numbers 
of spectator nucleons and on charged particle multiplicity 
are necessary to select collisions which produce sources 
with well-defined and large spatial deformation; the his- 
tograms shown in the bottom panel of Figure 21 corre- 
spond to only 0.025% of all U+U collisions taking place 
in the accelerator. The top panel in Figure 01 shows that 
it is possible to loosen these tight cuts somewhat, at the 
expense of reducing the average spatial source deforma- 
tion and introducing larger event-by-event fluctuations as 
well as an increased sensitivity to details of the Glauber 
model used for relating the relative nuclear orientation to 
the observed spectator nucleon and charged hadron mul- 
tiplicities. We leave a further discussion of such model 
uncertainties (see also |9j) for later when U+U collisions 
become (hopefully) available. 
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